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Abstract. The Bak-Sneppen model is an abstract representation of a biological system that 
evolves according to the Darwinian principles of random mutation and selection. The species 
in the system are characterized by a numerical fitness value between zero and one. We show that 
in the case of five species the steady-state fitness distribution can be obtained as a solution to a lin¬ 
ear differential equation of order five with hypergeometric coefficients. Similar representations for 
the asymptotic fitness distribution in larger systems may help pave the way towards a resolution of 
the question of whether or not, in the limit of infinitely many species, the fitness is asymptotically 
uniformly distributed on the interval {fc, 1] with fc ^ 2/3. 


E Introduction 

The Bak-Sneppen (B-S) model is an abstraet representation of a biologieal system that evolves 
aeeording to the Darwinian prineiples of random mutation and natural seleetion. It was intro¬ 
duced in [2] in the eontext of self-organized eritieality in systems with spatial interaetions. 

Despite its simplieity, the B-S model eaptures some of the features that are believed to be 
eharaeteristie of evolving biological systems. In particular, it predicts evolutionary activity on 
all time scales with long periods of relative stasis interrupted by bursts of activities, referred 
to as avalanches. As a eonsequenee of the absenee of a eharaeteristie time scale, evolutionary 
dynamies in the B-S model display long-range dependenee in both the temporal and the spatial 
domain. It is thus suitable as an abstraet representation of systems in punetuated equilibrium, a 
eoneept that was introduced in [6] to explain the patterns observed in fossil reeords. Eor a more 
thorough diseussion of these ideas, and for applieations of the notion of punctuated equilibrium 
in other seientifie disciplines we refer the reader to [12]. 

In addition to its usefulness in abstraetly representing some key features of palaeontology and 
maero-evolution, the B-S model has also been employed to analyze the evolution of baeteria in 
a eontrolled, eompetitive environment. In a series of key experiments [14], Lenski and eollabo- 
rators eultivated twelve initially identieal populations of an E. eoli strain over several years and 
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conserved samples at regular time intervals. They then determined the relative fitness of the eon- 
served samples by putting them into direet eompetition with a sample taken from the initial pop¬ 
ulations and measuring their relative growth rates. In [5] it was shown that the B-S model with 
random mutations qualitatively reproduces some of the experimental results on relative baeterial 
fitness obtained in Lenski’s long-term experimental evolution project [13]. Using an extended 
multi-trait variant, [3] extended the explanatory scope of the B-S model to inelude experimental 
findings about the interplay of adaptation, randomness and history in baeterial evolution. 


Informal definition of the B-S model. The Bak-Sneppen model eharaeterizes eaeh speeies in a 
biologieal system by a numerieal fitness value between zero and one, whieh represents its degree 
of adaptedness to its environment and ehanges as the speeies evolves. Further, eaeh speeies is 
assumed to direetly interaet with exactly two other species, where it is left unspeeified if such an 
interaction represents competition for resources, predator-prey relations or something different 
entirely. The Bak-Sneppen model can therefore be visualized as points on a eircle, where eaeh 
point stands for an eeological niche (or a species occupying that niche) and neighboring speeies 
interaet with eaeh other. 

As time progresses, the fitness landscape evolves in accordance with the following rules, rep¬ 
resenting in an abstraet way the principles of random mutation and natural seleetion: at eaeh 
time step, the least adapted species, i. e. the one with the smallest fitness parameter, is removed 
from the system (beeomes extinet) and its plaee is immediately taken by a new speeies whose 
fitness is initially modelled as an independent uniformly distributed random variable. In order to 
take into aeeount the effeet of this ehange on the loeal environment, the fitness parameters of the 
two speeies to either side of the least fit one are also reset to random values. This ean be thought 
of as those two speeies themselves beeoming extinet and superseded by new ones, or as them 
undergoing mutations in response to their neighbour beeoming extinet. 


Previous mathematical results. Despite its apparently easy definition, the B-S model has with¬ 
stood most attempts at mathematieal analysis in the past. Partial results have been obtained, 
however, in the eontext of rank-driven proeesses [10, 11] and mean-field approximations [4, 7]. 
Early on it was eonjeetured based on simulations that the steady-state fitness distribution at a 
fixed site eonverges, in the limit of large populations, to a uniform distribution on the interval 
[fc, 1], where fc is approximately equal to 0.667, but believed to be slightly larger than 2/3. 

There is only a small number of mathematieally rigorous result about the Bak-Sneppen model; 
in [15] it is shown that the steady-state fitness at a fixed site is bounded away from one in ex- 
peetation, independent of the number of speeies in the system; A eharaeterization of the limiting 
marginal fitness distribution, eonditional on a set of eritieal thresholds, is given in [16] (see also 
[8, 9]). In [18], the author proposes to eompute the steady-state fitness distribution as the fixed 
point of the one-step transition equation and uses this method to deseribe the asymptotie fitness 
distribution for four speeies in terms of a compaet rational function. In the same paper it is 
shown that one eannot find a similarly simple formula in the B-S model with five speeies, and 
that the fitness distribution of a randomly seleeted speeies at steady-state in this case is not only 
not rational, but not even a hypergeometric function. 
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Our contribution. In this paper, we revisit the Bak-Sneppen model with five speeies. In The¬ 
orem 1, our main theorem, we establish a representation of the steady-state fitness distribution 
for five species in terms of the solution of an explicit differential equation with hypergeomet¬ 
ric coefficients. This steady-state distribution encodes information about the fitness attributes of 
species in a system that has evolved for a long time. For instance, one can deduce from it how 
fit, on average, a randomly selected species from the population is expected to be; this is done in 
Corollary 1 . Furthermore, since the steady-state distribution contains information about the joint 
fitness values of all species in the population, its knowledge allows to draw biologically relevant 
conclusions about qualitative properties of the system, such as the emergence of one or several 
dominant species, or the fragmentation of the eco-system into areas of different prevailing fit¬ 
ness. In our simple model, the symmetry of the initial configuration is preserved and no such 
phenomenon occurs. It is an interesting question whether in systems with a more complicated 
interaction between species, symmetry can be spontaneously broken. 

We envisage that similar representations for the asymptotic fitness distribution in larger sys¬ 
tems may help pave the way towards a resolution of the question of whether or not, in the limit 
of infinitely many species, the fitness is indeed asymptotically uniformly distributed on the in¬ 
terval [/c, 1] with fc ^ 2/3. We speculate that the techniques developed in this paper can be 
generalized to analyze larger systems with more than five species. It seems plausible to expect 
that the asymptotic joint fitness distribution in such systems can still be characterized as the so¬ 
lution to a certain linear differential equation, even though the coefficients might no longer be 
hypergeometric functions. 


2. Formalization and main result 


We adopt the following formalization of the Bak-Sneppen model from [18]. Initially, all 
fitness parameters are independent uniformly distributed and after k evolutionary steps the state 
of the system is represented by the vector fj, 6 [0,1]^, where the fth component refers to the 
fitness of the zth species. The evolutionary dynamics of the system can be expressed formally by 
the equation 



where the one-step transition kernel Px encodes the dynamics of the model and is given by 






Here, and in the following, all vector indices are taken modulo five. The sequence f = (fk)k is a 
uniformly ergodic Markov chain with absolutely continuous marginal distributions with densities 
8k '■ [Oj 1]^ ^ This means that for any Borel set A G 1]^), 
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and that the random vectors converge in distribution to a steady-state limit foo. Moreover, the 
k-step densities gk satisfy the recursion 







( 1 ) 


where the vectors X]v[ 6 [0,1]^ and X]y[^ 6 [0,1]^ are obtained from x by dropping the vth, 
and (v + l)th components, or replacing these components by the components of respectively. 
Uniform ergodicity of the Markov chain f implies that the densities gk converge uniformly to the 
density g = goo of the unique invariant distribution of f, which we recognize as the steady-state 
fitness distribution. We also introduce the notation 


^n,;n(-T) — 2 ■ 




n,me 


X 6 


where 2 F 1 denotes the Gauss hypergeometric function [1, Section 15.1] and i = is the 
imaginary unit. The following is the main result of the paper. 


Theorem 1. The limiting density g = lim^_^oo gk given by 

5 

g(x) = l[oj]5(x)^<2'(min{Yy,Yy+i},max{Ay,Yy+i}), (2) 

V=1 


where q(x,y) = ^'(1 - x)S\{\ - y) + 5o,o(v). Here, 

Q{x) = ^F 2 J {1/2} Fi ,2 |yV2) + {1/2} F2,4 |yV2) , (3) 

and the function Si is the unique solution of the differential equation Y 4 ]=QCj{y)B^(^{y) = 0 with 
boundary conditions 

Si(l) = 1/5, s;(l) = 0, S"(l) = -1/5, sf (1) = 1, s;"\l) = -18/5. (4) 

The coefficients cfy), 7 = 0,1,..., 5, are hyper geometric functions given by 

Coiy) [y (y' - 22) &'(y) + (5y' - 14) Qiy)] , (5a) 

ci(y) =-ycoCy), (5b) 

C2(y) =-^ [y (3y' - 38y3 - 4) g'(y) + (l5y^ - lOy' + 4) ^(y)], (5c) 

C3(y) = - 12y [y (4y3 - l) ^'(y) + (Sy' + l) ^(y)], (5d) 

C4(y) = - 3y^ [y (y^ -h 2) Qfy) + (9y^ - 2) ^(y)], (5e) 

C5(y) =y' (y' + 2) [yQ'iy) - g(y)] . (5f) 

Finally, the function Bo a is given by 

1 

Bofl(x)= -[g'xome-Qiosfmdf. (6) 

Jl-x S 
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(a) Plot of the densities of the one-dimensional 
marginal distributions of for A: = 0 ,..., 6 
(dashed), together with their limit (solid line) 
as given by Equation (7). 



(b) Plot of the cumulative distribution func¬ 
tions of the marginal steady-state fitness dis¬ 
tribution in the Bak-Sneppen model with three 
(solid), four (dashed) and five (doffed) species, 
as well as fheir conjecfured limif as fhe number 
of species goes fo infinify (dash-doffed line). 


Unfortunately, it does not seem possible to evaluate the differential equation 2^=0 = 

0 in terms of known special functions. As a direct consequence of Theorem 1, we obtain the fit¬ 
ness distribution of a single species at steady-state by computing the one-dimensional marginal 
of Equation (2). 


Corollary 1. In the Bak-Sneppen model with five species, the fitness distribution of a randomly- 
selected species at steady-state is absolutely continuous with density 


Smarg(^^l — 




l[o,i]W- 


(7) 


Proof. Without loss of generality we can compute the distribution of the fitness of a fixed species, 
say the first one, which we call foo,\- Averaging over the fitness values of the remaining species 
and using Equation (2) leads to the expression 


wg(T) = r g{x, ^ 2 , ■ ■ •, = 1^ f 2 r 

j[0,lp \J Jo Jx 


l[0,i](x) 


for the density of /„ j. Plugging in the explicit formula for q stated in Theorem 1, differentiating 
once with respect to x and using Equation (6) to eliminate results in the equation g[y^^g(x) = 
Sj'(x). The observation that gmarg(O) = 3/5 together with Si(l) = 0 completes the proof. □ 


The graphs of the limiting density g^arg as well as the marginals of the k-step densities gk are 
depicted in Eigure la, illustrating the convergence asserted in Theorem 1. Moreover, Eigure lb 
compares the cumulative distribution functions of the steady-state fitness value at a fixed site in 
die Bak-Sneppen model with three (trivial), four [18, Theorem 1] and five (Theorem 1) species, 
as well as their conjectured limit - assuming /c = 2/3 - in a system with an infinite number of 
species. Clearly, a lot of work remains to be done to fill the gap between five and infinitely many 
species. 
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3. Proof of the main theorem 


It was shown in [18, Proposition 2] that the density of the joint distribution of the fitness 
parameters after k steps, starting from a uniform distribution 'l/([0,1])®^, is a polynomial given 
by 




■ 5 

^ qk (min {Xy, Xy+i], max {Xy, Xy+i }) 

-y=l 


where the funetions qk 6 Q[a, y] ean be written as 

qk{x,y) = ^ aij^kx'yk 

i,j>0 


The coefficients of q^ for k = 1,..., 5, as computed from Equation (1), are tabulated in ?? la-4. 
In order to identify the function g featuring in Theorem 1, it is thus sufficient to compute the 
uniform limit 

q(x,y) := lim qkix,y) = lim V aij^k^y. (8) 

k—^oo k—^oo f ^ 

ij 

This will occupy most of the rest of this section. First we recall an explicit recursion for the 
coefficients aij^k that was derived in [18, Proposition 5]. It corresponds to - and is derived from 
- the recursion (1) for the densities gk- It obviates the need to evaluate any integrals and thus 
allows for the functions gk to be determined much more quickly and efficiently. 


Proposition 1. The coefficients aij^k vanish for i = 0 and have the following properties: 
i) For j = 0, they satisfy aig^k = 0, (T 2 ,o Ml =2^1%' as well as 


CCifiMi - 1 + 

i + 2 


1 


+ 


rE 


O:i-2,0,k + 

A 
4 

/ Zj 


1 1 
3 ^ (i - 2)i 


ii - 1)/ 

. i JKTl . ^ 

I + 4 0^i-2,p,k I + 2 


Oli-3,0,k 


* n 1 2 / ^ + 1 

p=0 ^ p=0 ^ 


^i—2—p,p,k 


I ^ p + l 

p=0 ^ 


rZ 


2i ^ P 

p=0 ^ 

ii) For j > 1, they satisfy the recursion 

fySi+l o^kp.li . yj- 
I ^p=0 n+1 ' ^n=i 


, (-3 i-2 , i-3 

I + yp o:i-3-p,p,k yp cri-2-p,p,k 1 p _ 

i Zj p + 1 Zj i - p 2 Zj i - p " '' 


p=Q 


p=0 


0 P+1 + 2jp=o (p+\)(j-p)0:j-i-p,p,k^ 
(XiJ^k+l — ^IXijk 2®LL^+1’ 


i = 1, 

1 = 2, 


[o:i-ij,k - [l + + [I + a;iMj,k, i > 3. 


(9) 


( 10 ) 


The analysis of this three-dimensional recursion is simplified considerably by the fact that 
for each i, j the sequence (aij^k becomes eventually constant. More precisely, we have the 
following result which may be compared to [18, Lemma 2] and can be proved along the same 
lines. 
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Lemma 1. For each i, j there exists a rational number fii j such that = /3ijfor all k> i+j+ 1. 
In particular, lim^-^oo = Aj- 

In ?? la-4 the coefficients atj^k with k> i+j+ \ are printed in bold and the frontier coefficients 
aijj+j-i-i are marked with boxes. The assertion of Lemma 1 is thus easily appreciated by visual 
inspection. Equation (8) in combination with Lemma 1 allows to recognise the limiting function 
q as the generating function of the limiting coefficients /3ij, i. e. q(x,y) = Jjijfiijx'yj =: 5o,o(Y,y). 
Before proceeding further we define convenient notation for the generating functions of the ar¬ 
rays ajj^k and jSij along various dimensions and with various indices held fixed. 


Bojix) = ^ AjA, A,o(y) = ^fiijy\ 

i=0 j=0 

oo 

^i,j,o(.Z^ — CTjjyZ , ^o,j,ol_X,Z) — ^ Z , ^i,o,o(y^z) — ^ 0^i,j,ky^Z • 

k=0 i,k>0 j.k>Q 

Here, the formal variables x, y and z correspond to indices i, j and k, respectively, and the symbol 
o indicates summation over the index that it replaces. The next result establishes how passing to 
the limit k ^ oo can be accomplished at the level of generating functions. 

Lemma 2. For each x,y e [0,1], it holds that Bo^o(x,y) equals lim^^i- (1 - z)Ao,o,o(Y,y,z). Simi¬ 
larly, for non-negative integers i, j, it holds that 

Bi,o(y) = lim (1 - z)Ai^o,o{y,z), Bojix) = lim (1 - z)Aoj,o(x,z). (11) 

1 z^ 1 

Proof. It suffices to prove the first claim, which follows directly from Lemma 1 ; it allows us to 
write 

Ao,o,o(Y,y,z) = ^ ^ '^fiijizxYizyf. 

i,j k^i+j ij 

After multiplication by (1 - z) the first term vanishes as z approaches one, whereas the second 
one converges to Bo,o(Y,y). □ 


We now begin analyzing the recursion for or, in more detail. Throughout, we employ the 
powerful technique of generating functions as described in [19]. In particular, we make use of 
the fact that a linear recursion equation for a sequence c„ can be transformed into a differential 
equation for the generating function x by multiplying the original recursion by x” 

and summing over n. Indeed, if the original recursion equation has polynomial coefficients, this 
property is shared by the resulting differential equation. 


Proposition 2. For positive integers j, the generating function B^j is given by B^jix) = /3ijG(x), 
where 


GW = I 


i^4,5{l/2}F2,i 


(1 -xf 


-(1-y)2F2,i{1/2}L4,5 




( 12 ) 
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Proof. The last case of Equation (10) and Lemma 1 imply that, for j > I, the sequence 
satisfies the recursion 


fiij = Pi- 


L; 


1 + 


1 


i{i - 1) 


A-2,; + 


1 1 
3 ^ i{i - 2) 


> 3. 


After multiplying by iii - l)(z - 2)x' and summing over i this translates into the differential 
equation 

- 7(1 - x)B: .{x) + 3(1 - xfB'Pix) - ^ [2 + (1 - xf] Bf .^x) = 0 (13) 

for the generating functions B^jix) = Z”o The initial conditions are 


BojiO) = fioj = 0, B'^jiO) = /3,j, and 5"/O) = 2/32j = Pi 


j’ 


(14) 


where the last equality follows from the second case of Equation (10). The general solution of 
Equation (13) is 


Bojix) =Cojix - 1 ) 3 E 2 (l, i (3 - i V 2 ), ^ (3 + i V 2 ); ^; ^(Y - 1 ) 


+ Clj(l “ XpF4 5 


(l-xf 


+ C2jF2,3 


3 3 2 
(1-y)3 


This can be obtained with the aid of a computer algebra system or checked using the power series 
representation of hypergeometric functions. The coefficients c,nj, m = 1,2,3, are determined by 
the initial conditions (14) and are given by cqj = 0, Ci,^ = diPij, and C 2 ,j = d 2 P\,j, where d\ 
and d 2 are explicit hypergeometric constants. Observing that di + d 2 equals 40/9, which can be 
proved via Zeilberger’s algorithm [17], the claim follows. □ 


Proposition 2 implies the decomposition 

00 

Bo,o(x,y) = ^ Boj(x)y^ + 5o,o(y) = G(x)Bio(y) + B^flix), 

;=i 

and it thus only remains to compute the two functions Bi o and Bofl. Unfortunately, the corre¬ 
sponding cases of Equations (9) and (10) can not be analyzed by simply passing to the limit 
k ^ oc because they involve coefficients a, with j > k. We therefore need to compute the 
two-dimensional generating function of the array a task which is directly modelled after 

the proof of Proposition 2. 


Proposition 3. For positive integers j, the generating function Aoj^o'- {x,z) is 

given by Aoj^o{x,z) = Ai jo(z)U 2 (y,z), where 


G2{x, z) 


9 

2(3-zf 



3-z / 



(1 

3-z j ■ 


( 15 ) 


In particular, G 2 (x, 1) = G(x). 
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Proof. The proof proceeds by multiplying the i > 3 case of Equation (10) by i{i - !)(/ - 2)x'z^ 
and summing over both i and k. After some algebraic manipulations one arrives at a differential 
equation similar to Equation (13) which can be solved using the same techniques. □ 

Our next task will be to transform the first case of the recursion relation (10) into differential 
form. This corresponds to the case y > 1 and z = 1. The proof proceeds as before by multiplying 
the recursion equation by summing over j and k and interchanging the order of summation 
where necessary. 

Proposition 4. The generating functions Ai o,o and Ao,o,o satisfy the integral equation 

- (^i,o,o(y,z) - 1) =Giy,z) f (Ai,o,o(vp,z) - I)!)!!- + (Ai,o,o(y,z) - 1) 

Z Jy 

+ (1-y)Ao,o,oCv,z) + r Ao,o,o(w,z)dw. 

Jo 

Setting z = 1 in Equation (16), applying Eemma 2 and differentiating once with respect to x 
yields the following result which will be one of the main ingredients in the proof of our main 
theorem. 


f 


G(w, z)dvp 


(16) 


Corollary 2. The generating functions Bi o and Bog satisfy the integro-dijferential equation 

>1 


BUx) 


1 


\ - X 


-G'(x)£ BUf)d^ + B[^(x) h: G(f)df 


(17) 


A second differential equation relating Bi o and Bog can be derived from Equation (9). 


Proposition 5. The generating functions Bgo and Bog satisfy the differential equation 
4(1 -x)£ Bog(f)df - (1 - x)^Bog(x) + 1 [2 + (1 - y)3] B'fyx) 

= 3(1 - ^)(l - G(f)dfjBUx) + ^ [2 + (1 - X)'] G(x)B' o(x). (18) 

Proof Using the i = 1 case of Equation (10) the two ZptV sums can be eliminated from Equa¬ 
tion (9). This allows to pass to the limit k ^ oo and to replace the coefficients atj^k by their 
steady-state values fiij from Eemma 1. Multiplying the resulting equation by i(i - l)(z - 2)x‘ 
and summing over i yields Equation (18). The details of the computation are lengthy and hence 
omitted. □ 


The two equations (17) and (18), together with the initial conditions that are tabulated in 
?? la-4, are already enough to numerically compute the functions Bi o and Bog. It is, however, 
possible to simplify the problem even further to a single differential equation of order five. 


Proof of Theorem 1. Eor better readability we introduce the antiderivatives 


So(y) 


r 


Bogimf, Si(x) 


i: 


Bifyf)df, and ^(x) = 1 


-f' 


G(Odf. 
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A direct evaluation of the last antiderivative gives Equation (3). After the substitution x ^ ^-y. 
Equations (17) and (18) can be rewritten in terms of these functions as 

yK(y) = - ^''(y)Si(y) + Q(y)!B';(y), 

4y!Bo(y) + y^S'^iy) + ^ (2 + Sg (y) =3yQ(y)!B[(y) - ^ (2 + Q'(y)S'((y). ^ ^ 

Differentiating the latter equation three times eliminates the zeroeth and first derivative of So; 
replacing the remaining second through fifth derivatives with expressions derived from the first 
equation of the last display produces the fifth-order differential equation Ci(y)S^f(y) = 0. 
The coefficients cj are given by 

co(y) =6 (5y' - 2) g''(y) + 6y (Sy^ + 2) + y^ (9y^ - b) 0^%) + y^ (y^ + 2) 

Ciiy) =3y [(y' + 4) ^"(y) + y (3y' - 4) + y^ (y^ + 2) ^(^^(y)], 

C 2 (y) =6 (2 - 5y') ^(y) - 6y (l3y' + 2) g'(y) - 9y^g"(y) + y' [(l ly' + 4) g<-\y) + (y^ + 2) yg^%)], 
csiy) =y [-3 (l9y' + 4) g(y) + 3y (4 - 9y') ^'(y) + 4y^ (4y^ - l) g"{y) + 3y' (y' + 2) g^^^(y)], 

C 4 (y) =3y^ [(2 - 6y^) 0{y) + 2y (y^ - l) g'(y) + y^ (y^ -P 2) 0"(y )], 
c5 (3^) =3^^ + 2) (-^(y) + y^'(y)). 

It follows from Proposition 2 that the function ^ satisfies the differential equations 

0 =3yg(y) + 3yV(y) + (y' + 2)g"(y), 

0 =6y^0'(y) + (5y' - 2) g"(y) + y (y' + 2) ^('^(y), 

0 = (l lyH 4) g"(y) + y (ly^ - 4) g^\y) + y^ (y^ + 2) ^("^(y), 

0 =18y (l ly^ + 16) g^^'>(y) + 9y^ (l ly^ - 2) ^(^^(y) + (l ly^ + 26y3 + 8) ^(^^(y). 

Using these equations to eliminate the higher-order derivatives of ^ from the expressions for 
the coefficients Cj(x) one obtains the claimed expressions 5. Assertion (6) is a reformulation of 
Equation (17). □ 
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Appendix A. Coefficients 

In this appendix we eolleet the values of the eoeffieients a^j^k of the funetions qt for k = 
1,..., 5. This serves as an illustration of Lemma 1 and to determine the initial eonditions of the 
various differential equations eneountered in the paper. 


\j 

i \ 

0 1 

0 

0 0 

1 

1 0 

2 

-1 0 

3 

1/3 0 

4 

0 0 


{A)k = \ 



(b) k - 2 


Table 1. Coeffieients of gu, k = 1,2, as eomputed from Equation (1) or 
Equations (9) and (10). 



0 12 3 


0 

0 


0 

0 

0 

1 

0 


1/5 

29/10 

-25/3 

2 


2/5 


9/10 

-59/20 

31/6 

3 

38/5 


-5/3 

5/2 

-5/3 

4 

-523/20 

1 

-3/2 

1 

5 

75/2 


-1/5 

3/10 

-1/5 

6 

-3551/120 

0 

0 

0 

7 

477/35 

0 



8 

-487/140 

0 



9 

487/1260 

0 



10 

0 


0 




4 

5 

6 

7 

8 

0 

0 

0 

0 

0 

121/12 

-32/5 

32/15 

-32/105 

0 

-127/24 

16/5 

-16/15 

16/105 

0 

5/12 

0 

0 

0 

0 

-1/4 

0 




1/20 

0 




0 

0 





Table 2. Coefficients Qr, j 3 of ^3 as computed from Equation (1) or Equations (9) 
and ( 10 ). 
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j 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

I 

0 

1/5 

03 

7 

-1879/60 

8507/150 

-10421/180 

4589/126 

-227/16 

143/45 

-143/450 

0 

2 

2/S 

00 

53/20 

-71/6 

3089/120 

-10427/300 

11189/360 

-23329/1260 

227/32 

-143/90 

143/900 

0 

3 

0 

2/3 

-19/3 

134/9 

-307/18 

32/3 

-32/9 

32/63 

0 

0 

0 

0 

4 

1019/60 

-51/20 

281/40 

-133/12 

517/48 

-32/5 

32/15 

-32/105 

0 

0 

0 

0 

5 

-2591/30 

311/100 

-1061/200 

289/60 

-673/240 

32/25 

-32/75 

32/525 

0 

0 

0 

0 

6 

39877/225 

-223/120 

223/80 

-223/120 

223/480 

0 

0 

0 

0 

0 

0 

0 

7 

-138271/2100 

17/30 

-17/20 

17/30 

-17/120 

0 

0 

0 

0 

0 

0 

0 

8 

15035/96 

-17/240 

17/160 

-17/240 

17/960 

0 

0 

0 

0 

0 

0 

0 

9 

-55459/720 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 

1224179/50400 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

11 

-113287/25200 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

12 

113287/302400 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Table 3. Coefficients or, y _4 of g 4 as computed from Equation (1) or Equations (9) 
and (10). 
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0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1/5 

1/2 

Ezg 

319/20 

-4889/50 

71923/300 

-357143/1050 

879181/2800 

-62029/315 

53009/630 

-12424/525 

896/225 

-896/2925 

0 


2/5 

1/10 

EEi] 

67/10 

-943/24 

31681/300 

-319979/1800 

1300879/6300 

-958631/5600 

64031/630 

-267047/6300 

6212/525 

-448/225 

448/2925 

0 


1 

ElE 

31/15 

-20 

1121/18 

-9083/90 

53257/540 

-115301/1890 

1135/48 

-143/27 

143/270 

0 

0 

0 

0 

4 

147/601 

13/20 

-359/40 

371/12 

-14231/240 

11147/150 

-11477/180 

23473/630 

-227/16 

143/45 

-143/450 

0 

0 

0 

0 


5651/150 

-321/100 

2947/200 

-629/20 

46771/1200 

-23627/750 

16469/900 

-25969/3150 

227/80 

-143/225 

143/2250 

0 

0 

0 
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-37733/150 
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-17927/1200 
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-5855/288 
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892/1575 
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0 

0 

0 
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0 
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36457/5040 
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0 

0 

0 

0 

0 

0 
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9967/2800 

-93511/16800 
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-32411/20160 

34/75 
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34/1575 

0 

0 

0 

0 

0 

0 
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-1957/1512 

1957/1008 

-1957/1512 

1957/6048 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 

-20505517/25200 

289/1080 

-289/720 

289/1080 

-289/4320 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

11 

694474463/1663200 

-289/11880 

289/7920 

-289/11880 

289/47520 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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^97946013/3326400 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

14 
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0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

15 

6257393/18018000 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


Table 4. Coefficients of gs as computed from Equation (1) or Equations (9) 
and (10). 










